Heterogeneous spin state in the field-induced phase of volborthite as seen via ^^V 

nuclear magnetic resonance 
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We report results of '^^V NMR in the field-induced phase of volborthite Cu3V207(OH)2-2H20, a 
spin-1/2 antiferromagnet on a distorted kagome lattice. In magnetic fields above 4.5 T, two types 
of V sites with different spin-echo decay rates are observed. The hyperfine field at the fast decaying 
sites has a distribution, while it is more homogeneous at the slowly decaying sites. Our results 
indicate a heterogeneous state consisting of two spatially alternating Cu spin systems, one of which 
exhibits anomalous spin fluctuations contrasting with the other showing a conventional static order. 

PACS numbers: 75.30.Kz, 76.60.-k, 75.40.Gb 



Possibility of exotic quantum states in two-dimensional 
(2D) spin systems with frustrated interactions have at- 
tracted strong attentioniii^ In particular, the ground 
state of the spin-1/2 Heisenberg model with a nearest- 
neighbor antiferromagnetic (AF) interaction on a kagome 
lattice, a 2D network of corner-sharing equilateral tri- 
angles, is believed to develop no long-range magnetic 
order. Theories have proposed various ground states 
such as spin liquids with no broken symmetry with^i 
or without^ an excitation gap or symmetry breaking 
valence-bond-crystal statesi^ Candidate materials known 
to date, however, depart from the ideal kagome model 
in some aspects, such as disorder, structural distor- 
tion, anisotropy, or longer range interactions. Effects of 
the Dzyaloshinsky-Moriya (DM) interaction,— spatially 
anisotropic exchange interactions ji"— and longer range 
Heisenberg interactionsiS have been theoretically inves- 
tigated. 

Volborthite Cu3V207(OH)2-2H20 is a layered com- 
pound, in which distorted kagome layers are formed by 
isosceles triangles of Cu ions carrying a spin 1/2. Con- 
sequently, there are two Cu sites, Cul and Cu2, and 
two kinds of exchange interactions, J and J', as shown 
in Fig. 1. The Cu2 sites form linear chains, which are 
connected through the Cul sites. The magnetic suscep- 
tibility X obeys the Curie- Weiss law x — C/{T + 6w) 
above 200 K with 9w = 115 K, exhibits a broad maxi- 
mum at 20 K, and approaches a finite value at the lowest 



temperatures ji 



An unusual magnetic transition has 



been observed near 1 K^^— which is much lower than 
6w, consistent with strong frustration in a kagome lat- 
tice. However, a recent density functional calculatioiiii 
proposed a ferromagnetic J and a sizable AF interaction 
J" between second neighbors along the chain as shown 
in Fig. 1. In this model, frustration arises from the com- 
petition between J and J" along the chain rather than 
the geometry of a kagome lattice. The appropriate spin 
model for volborthite has not been settled yet. Recently, 
anomalous sequential magnetization steps were reported 
in a high quality polycrystalline sample at 4.3, 25.5, and 
46 TJ^ Kaneko et al. showed a classical Heisenberg 



model on a spatially distorted kagome lattice exhibits 
a magnetization step by the Monte Carlo methodJ^ 

A magnetic transition in volborthite at T* ~ 1 K 
was detected by ^^V-NMR )^^i^^ muon spin relaxation^i^ 
and heat capacityi^ experiments. Our previous NMR 
measurements^ revealed a sharp peak in the nuclear 
spin-lattice relaxation rate 1 /Ti and onset of broadening 
of the NMR spectrum due to development of spontaneous 
moments for the magnetic field B below 4.5 T. A kink was 
also observed in the heat capacity^i^ Most recently, devel- 
opment of short range spin correlation was detected by 
neutron inelastic scatteringJ^ However, the low T phase, 
which we call phase I, shows various anomalies incompat- 
ible with a conventional magnetic order The NMR 
line shape is not rectangular but fits to a Lorentzian 
(Fig. 2b), suggesting spatial modulation of static mo- 
ments. The behavior 1/Ti oc T provides evidence for 
dense low-energy excitations. The spin-echo decay rate 
l/r2 is anomalously large, pointing to unusually slow 
spin fluctuations. 

Above 4.5 T, at which the first magnetization step oc- 






FIG. 1: (color online) Schematic structure of volborthite pro- 
jected onto the a-b plane. The H and O sites are not shown. 
The V sites are located below and above the Cu kagome lay- 
ers. The upper and lower V sites are related by inversion with 
respect to the Cu sites, hence all V sites are equivalent. See 
the text for the definition of various symbols. 
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curs, another magnetic phase appears with the different 
NMR line shape and dynamicsJ^ The spin state in tliis 
phase, which we caU phase II, has not been well charac- 
terized yet. In this paper, we report results of ^^V NMR 
in phase II, which revealed coexistence of two types of V 
sites with different I/T2 and line shapes. Our results in- 
dicate that phase II exhibits a unique heterogeneous spin 
state, where two distinct types of Cu spins, one with non- 
uniform magnitude and large temporal fluctuations and 
the other with a more conventional static order, form a 
periodic structure. 

^^V NMR measurements have been performed on a 
high-quality powder sample synthesized as described in 
Ref. The NMR spectra were obtained by summing 
the Fourier transform of the spin-echo signal at equally 
spaced frequencies in a fixed B. The pulse sequence 
tt/2 — T — 7r/2 was used with the pulse width of 1-3 
/US and the pulse separation r in the range of 10-150 
IIS. We determined I/T2 by fitting the spin-echo in- 
tensity /(r) as a function of r to an exponential decay 
/(r) = AoCxp{—2t/T2). Likewise, 1/Ti was obtained by 
fitting the nuclear recovery curve /(i), which is the spin- 
echo intensity as a function of the time t between the 
saturating comb pulses and the first 7r/2 pulse, to the 
exponential function I{t) = leg — /oCxp(— t/Ti). When 
this function did not fit the data due to inhomogeneous 
distribution of l/Ti, we used the stretched exponential 
function I{t) = I^q — /ocxp{ — (t/Ti)'^} to determine the 
representative value of 1/Ti. 

Figure 2(a) shows the r dependence of the NMR spec- 
tra at 6.0 T and 0.3 K (phase II). The spectrum for r = 

10 IjLB is identical to one of the spectra shown in Fig. 4(a) 
of Ref. [1^ except that it was plotted against the shift in 
resonance field rather than the frequency in Ref. [l^ In 
addition to the overall decay of the intensity with r, the 
line shape also changes with r. The broad peak near the 
center of the spectrum observed for r = 10 and 25 /xs dis- 
appears for r > 40 /xs. Instead the two peaks indicated by 
the arrows become more pronounced with increasing r. 
Such T dependence indicates that the spectrum consists 
of two or more components characterized by different val- 
ues of 1/T2. On the other hand, the line shape at 1 T 
(phase I) does not depend on r as shown in Fig. 2(b). 

We confirmed the two-component behavior in phase 

11 by direct measurements of I/T2. Figure 3 shows the 
spin-echo decay curves at 0.3 K measured at the center of 
the spectrum obtained in the field of 1, 6, and 12 T. The 
data at 1 T can be fit to the single exponential function, 
providing unique I/T2. On the other hand, the data at 
6 and 12 T show a convex curve, which can be well fit to 
the two-component exponential function 

/(r) = A/exp(-2r/r2/) + A,exp(-2r/r2,) (1) 

as shown by the solid lines. The fast and slowly decaying 
components at 6 T are shown separately by the dashed 
lines. The decay rate of the fast component, I/T2/; is 
nearly equal to I/T2 at 1 T. 
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FIG. 2: (color online) NMR spectra for various values of r 
at 0.3 K in the magnetic field of (a) 6 T and (b) 1.01 T. 
The spectra in (c) and (d) show decomposition of the NMR 
spectra at 0.3 K into the fast decaying (solid circles) and 
slowly decaying (open circles) components for the magnetic 
field of 6 T and 12 T. 



We define the two V sites showing the fast and slow de- 
cay rates as the V/ and Vs sites, respectively. The coeffi- 
cients Af and in eq. (1) are proportional to the num- 
ber of these sites within the frequency window covered by 
the exciting rf-pulse. Therefore, the decomposition of a 
NMR spectrum into the fast and slow components can be 
accomplished by measuring the spin-echo decay curve at 
many frequency points covering the entire spectrum and 
plotting Af and As against frequency. The decomposed 
spectra thus obtained at 6 and 12 T are shown in Fig. 2(c) 
and (d) . The pulse width was set to 3 ^s covering about 
±0.1 MHz. The two-component fit of the spin-echo decay 
curve was successful at all frequencies. The integrated in- 
tensities of the "fast" and "slow" spectra correspond to 
the total number of the V/ and Vs sites. For the spectra 
shown in Fig. 2(c) and (d), the integrated intensities of 
the two components are nearly equal within the experi- 
mental error, indicating equal abundance of the V/ and 
Ys sites. 

The inset of Fig. 3 shows the B dependence of I/T2 
at 0.3 K measured at the center of the spectrum. In 
phase I, the spin-echo decay curve can be always fit to 
the single exponential function. Such a clean behavior 
is brought by the improvement of sample quality)^ since 
two-component behavior was observed even near 2 T in 
earlier studies The decrease of I/T2 with increasing 
field indicates suppression of slow spin fluctuations by 
magnetic fields. On the other hand, 1/Ti is insensitive 
to B in phase IJ^ In phase II, the decay curve can be 
well fit to the two-component function (1). Both I/T2/ 
and 1/72.5 slightly decreases with increasing B. Near the 
transition at 4.5 T, the two-component fit docs not work 
well probably due to coexistence of the two phases. 
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FIG. 3: (color online) Spin-echo decay curves at 0.3 K in the 
magnetic field of 1, 6, and 12 T. The fast and slowly decaying 
components for the 6 T data are shown by the dashed lines. 
The inset shows the B dependence of I/T2 at 0.3 K. The open 
squares indicate I/T2 in phase I. The solid and open circles 
represent 1/T2f and I/T2S in phase II. 



Figure 4(a) shows the T dependence of I/T2 at 6 T 
together with the previous data at 1 TJ^ The spin-echo 
decay at 1 T follows the single exponential function at all 
temperatures (crosses). At 6 T, the decay curve can be 
fit by the single exponential function for 1.5 < T < 2 K. 
When I/T2 becomes small for T > 2.0 K, the Gaussian 
contribution to the spin-echo decay due to the nuclear 
spin-spin coupling has to be taken into account. Then 
the decay curve can be fit to the "exponential -I- Gaus- 
sian" function /(r) = Aoexp{-2T/r2 - 2(r/r2G)^}. The 
Gaussian decay rate is almost temperature independent 
(l/r2G = 2 xlO^ s~i). The decay curve exhibits the 
two-component exponential behavior below 1.5 K, which 
is the boundary between the uniform paramagnetic phase 
and phase IL We then determined 1/72/ (solid circles) 
and I/T2S (open circles) by using eq. (1). Although I/T2/ 
and I/T2S are different in magnitude, they show similar 
T dependences and decrease slightly with decreasing T. 
This T-dependence is nearly identical to the behavior at 
1 T (Fig. 4a). Below 1.1 K, the ratio AJAj is almost T- 
independent, indicating that the numbers of the V/ and 
Vg sites stays nearly equal. The two-component fit docs 
not work well near 1.5 K probably due to coexistence of 
the paramagnetic phase and phase IL 

Figure 4(b) shows the T dependence of 1/Ti at 6 T 
(triangles) together with the previous data at 1 T 
(asterisks) The recovery curve at 1 T follows the ex- 
ponential function for T > 1 K, while the stretched ex- 
ponential function is required to fit the data at lower 
temperatures^^ Although the recovery curves do not de- 
pend on T at 1 T (phase I) , they do in phase IL The open 
and solid triangles represent 1/Ti at 6 T for r = 10 and 
80 IIS, respectively. Above 1.5 K, the recovery curve can 
be reproduced by the exponential function and 1 /Ti does 
not depend on r. Below 1.5 K (phase II) the recovery 
curve for r = 10 /xs was fit to the stretched exponential 
function. We then obtained 1/Ti oc T^-^ behavior but 
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FIG. 4: (color online) The T dependences of (a) I/T2 and 
(b) 1/Ti. The various symbols are defined in the text. Two 
components of I/T2 appear at 6 T below 1.5 K, whereas only 
the slow component can be determined for 1/Ti. The data at 
1 T are taken from Ref. [Tsl . 



/3 depends on T and reaches an extremely small value of 
0.25 at 0.14 K. This can be now accounted for by coex- 
istence of the V/ and sites for short r. The recovery 
curve for r = 80 ^s can be fit to the stretched exponen- 
tial function with constant /3 « 0.6, which is also the case 
at 1 T. Since the signal from the V/ sites has decayed 
at long T and is nearly absent at r = 80 fj,s, 1/Ti thus 
obtained should correspond to the Vs sites. The solid 
line in Fig. 4(b) shows a 1/Ti (x T behavior, which fits 
the data for both 6 T and 1 T. This indicates that the 
Vs sites in phase II are affected by anomalously dense 
low-ener gy e xcitations that have been identified in phase 
I in Ref. [l^. The smaller values of l/Ti at 6 T indicates 
that the Vs sites couple to these excitations more weakly 
than the V sites in phase I. 

How can we understand the spin state in phase II? The 
V/ and Vs sites show significantly different relaxation 
rates and line shapes. The spectra at the V/ sites show 
a round shape (Fig. 2cd), indicating a distribution in the 
magnitude of the hyperfine field, which can be generated, 
for instance, by a spin-density-wave order or spatially 
disordered static moments. Therefore, the large 1/12/ at 
the V/ sites should be ascribed to unusually slow fluctu- 
ations of the Cu spins with non- uniform moments. On 
the other hand, the Vs sites show a rectangular spectral 
shape, which is compatible with conventional AF order. 
The relatively small l/r2s and 1/Ti show that the effects 
of the unusual spin fluctuations are weak at the Vs sites. 

How are the two V sites distributed spatially? Two 
cases can be considered. In one case, a magnetic super- 
structure is formed in phase II yielding two inequivalent 
V sites. Alternatively, macroscopic phase separation may 
occur and each phase contains only one of the two V sites. 
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Our NMR results definitely support the first case. The 
numbers of the and Yg sites are nearly equal in the en- 
tire region of phase II, which covers a wide range of B and 
T. This is quite natural if a microscopic superstructure 
is formed but very unlikely for the phase separation sce- 
nario, since the relative stability of the two phases should 
change with B and T. In addition, the two V sites show 
similar T dependences of I/T2 in spite of the difference 
in magnitude (Fig. 4a). This indicates a common source 
of the fluctuating hyperfine field for the two V sites. The 
different magnitude should then be attributed to the dif- 
ference in the hyperfine coupling between the fluctuating 
spins and the nuclei. Since the range of hyperfine inter- 
action is short with an atomic scale, this is not possible 
for the case of phase separation. 

Although we cannot determine the superstructure from 
the NMR data on a powder sample alone, let us discuss 
most likely possibilities. Since the numbers of the two V 
sites are equal, we expect a doubled unit cell due to for- 
mation of a superstructure. For the ideal kagome lattice, 
the three directions ki, k2, and ki -|- k2 (Fig. 1) are all 
equivalent. However, since the structure of volborthite 
is uniform along k2 and incquivalcnt Cul and Cu2 sites 
alternate along ki, we consider the wave vector of the 
superstructure is likely to be along ki. The two V sites 
then alternate along ki. The Vq, and in Fig. 1 cor- 
respond to the V/ and Vg sites. Likewise, each of the 
Cul and Cu2 sites are divided into two inequivalent sites 
Gulp. CuIct and Cu2k, Cu2\, as shown in Fig. 1. 

Since a V site is located approximately above or below 
the center of the Cu hexagon of the kagome lattice, the 
dominant source of the hyperfine field at V nuclei should 
be confined within the six Cu spins on a hexagon. Be- 
cause of the distorted structure of volborthite, there are 
three distinct hyperfine coupling tensors Aq, Af,, and 
Ac as shown in Fig. 1. The coupling tensors to the 
other three spins A'^, A'^^, and A^ are generated by the 
mirror reflection perpendicular to the b axis at the V 
site. The hyperflne flelds at Vq and Yp arc written as 
B„ = {AaS^+A'^s'^) + {AbSp+A'^s'p) + {AcSx+Ks'^) and 
Bfi = {AaSx+Ks\) + {AbS^+A'^s'^) + {A^s^+A'^s'^), re- 
spectively, where Se and denote two neighboring spins 
on the same type of sites, (e stands for Gulp, Culo-, 



Cu2„, or Cu2a.) 

Both 1/Ti and I/T2 are determined by the time- 
correlation function of the hyperfine field<^ If all the 
hyperfine coupling tensors A^, Af,, and Ac are largely 
isotropic and have similar magnitudes, the different l/r2 
for the two V sites must be ascribed to the difference in 
the time-correlation functions of and Sa-. A plausi- 
ble situation is illustrated in Fig. 1. The Culo- sites de- 
velop a conventional long range magnetic order. Hence, 
the spin fluctuations of these sites have a small ampli- 
tude with a conventional time scale typically described 
by spin waves. On the other hand, the Gulp sites show 
modulated spin structure. The absence of fully developed 
static moments would allow unusually slow fluctuations 
responsible for the large I/T2 at the V/ sites. By study- 
ing the anisotropic kagome model in the limit J ^ J', 
Schnyder et at proposed a spiral order with large (small) 
ordered moments on the Cul (Cu2) sitesii Although the 
different magnitudes of the ordered moments for the two 
Cu sites is a feature in common to our results, our obser- 
vation of two distinct V sites requires an enlarged mag- 
netic unit cell. The spin structure of the Culo- sites could 
be a spiral order with a small wave vector as proposed 
theoreticallyjii^ 

Alternatively, the low symmetry of the volborthite 
structure may result in large difference in the hyperfine 
coupling. In particular, A^ and Ac involve significantly 
different V-O-Cu hybridization paths. If one of them, say 
Aa, is the dominant coupling, the different I/T2 for the 
two V sites must be ascribed to the different dynamics of 
s\ and s„ spins. 

In conclusion, our NMR results show coexistence of 
the two types of V sites with different relaxation behav- 
ior and spectral shapes in the field-induced phase II of 
volborthite. Our results indicate that the high field phase 
of volborthite exhibits a unique heterogeneous magnetic 
state, where a non-uniform order with anomalous fluctu- 
ations alternates with a more conventional static order. 
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